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ABSTRACT 

The mean intracluster gas fraction of X-ray clusters within their hydrostatic re¬ 
gions is derived from recent observational compilations of David, Jones & Forman 
and White & Fabian. At radii encompassing a mean density 500 times the critical 
value, the individual sample bi-weight means are moderately (2.4cr) discrepant; re¬ 
vising binding masses with a virial relation calibrated by numerical simulations re¬ 
moves the discrepancy and results in a combined sample mean and standard error 
fgas{r 50 o) = (0.060 ± 0.003) For hierarchical clustering models with an extreme 

physical assumption to maximize cluster gas content, this value constrains the univer¬ 
sal ratio of total, clustered to baryonic mass < 23.1combining with the 

primordial nucleosynthesis upper limit on U;, results in < 0.60. A less conserva¬ 

tive, physically plausible approach based on low D/H inferences from quasar absorption 
spectra and accounting for baryons within cluster galaxies yields an estimate 

^ 0.28 ±0.07 

with sources of systematic error involved in the derivation providing approximately 
35% uncertainty. Additional effects which could provide consistency with the Einstein- 
deSitter case = 1 are presented, and their observable signatures discussed. 
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1 INTRODUCTION 

Clusters of galaxies provide a number of interesting cos¬ 
mological diagnostics. In particular, the relative amount of 
baryons and dark matter within their hydrostatic regions 
provides a measure of the cosmic mix of these components, 
i.e., a measure of the ratio of density parameters Qrn/^b iu 
the Friedmann-Lemaitre world modeE 

Employing this measure in practice requires accurate 
observational data along with estimates of possible system¬ 
atic biases. Likely sources of bias include systematic errors in 
component mass estimates and deviation of the local cluster 
ratio of baryonic-to-total mass arising, for example, from 
the different dynamical histories of the two components. 
Within the context of hierarchical clustering scenarios, these 
effects have been calibrated by numerical simulations of clus- 

* In this paper, Qrn refers to the contribution of all clustered 
matter (including baryons) to the stress-energy density. 


ter formation. The results, discussed in detail below, indicate 
that the magnitude of these effects are small — a few tens 
of percent or less — in observationally accessible regions 
of clusters. If our current description of cluster formation 
dynamics is physically accurate, then determination of the 
cosmic baryon fraction from X-ray cluster data is straight¬ 
forward. 

The baryonic component of the richest clusters is dom¬ 
inated by the X-ray emitting intracluster gas rather than 
the mass associated with the optical light of the cluster 
galaxies (Forman & Jones 1982; Sarazin 1986; Mushotzky 
1994). Taking Coma as an example. White et al. (1993) es¬ 
timate the ratio of gas to galaxy mass to be Mgas/Mgai = 
(5.5 ± 1.5) (with h = Ho/100 km s“^ Mpc”*^). This 

number is typical of the values seen in larger samples for 
clusters comparable in temperature to Coma (David et al. 
1990; Arnaud et al. 1992). The data also indicate a depen¬ 
dence of the ratio Mgaa/Mgai on cluster temperature, with 
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poorer clusters and groups possessing less gas per massive 
galaxy than their larger counterparts. 

The gas mass thus provides not only a formal lower limit 
to the total baryon cluster content, but a fairly accurate es¬ 
timate of the baryon mass in rich clusters if the dark mat¬ 
ter is assumed non-baryonic. Recent compilations of X-ray 
cluster data by White & Fabian (1995) and David, Jones & 
Forman (1995) provide gas and total mass estimates for sam¬ 
ples of moderate size. In this paper, I use these data, along 
with guidance from numerical simulations, to estimate the 
sample mean gas mass fraction fgas within a characteristic 
radius defining the hydrostatic boundary of clusters. Moti¬ 
vations for the specific approach adopted here are provided 
in §2 and the observational data are analysed in §3. Implied 
constraints on the universal baryon fraction and the value 
of Dm are provided in §4, including a thorough discussion of 
systematic errors. 


2 PHYSICAL EXPECTATIONS 

The value of any individual cluster’s baryon-to-total mass 
inferred from observations will, in general, differ from the 
cosmic value. The difference can be partly intrinsic — re¬ 
flecting a true baryon enhancement or deficit within the 
cluster — and partly due to errors in estimates of the com¬ 
ponent masses. 

To the extent that the physical processes responsible 
for a cluster’s structure are independent of its total mass, 
any intrinsic bias in the baryon fraction can, to first approx¬ 
imation, be expressed as a function of a scaled radial vari¬ 
able, equivalent to the local density contrast. Introducing 
notation used below. Sc denotes the mean interior density 
contrast with respect to the critical value Sc = p{< r)/pc 
with Pc = 3 Ho^/SttG. Similarly, rs^ is the radius at which 
a density contrast Sc is attained. 

2.1 Intrinsic component bias 

Because the horizon mass scale at the baryogenesis epoch 
is many orders of magnitude smaller than the typical clus¬ 
ter mass (e.g., Kolb & Turner 1990), there is no causal 
mechanism which can generate primordial fluctuations in 
the baryon-to-total mass ratio on cluster scales. This argu¬ 
ment holds as long as inflation precedes baryogenesis. Any 
intrinsic bias must therefore be set by dynamical processes 
operating differentially on the baryonic and dark matter. 

Gravity alone is not a powerful segregating mechanism. 
In the simplest example, consider clusters forming from 
spherically symmetric, scale-free initial density perturba¬ 
tions. Self-similar solutions to the dynamical equations (Fill¬ 
more & Goldreich 1984; Bertschinger 1985; Ghieze, Teyssier 
& Alimi 1996) exhibit two characteristic regions — a nearly 
hydrostatic, inner body surrounded by an outer, infalling 
envelope which merges seamlessly into the Hubble flow at 
large radii. The flow changes discontinuously at the bound¬ 
ary, implying the development of a shock for the collisional 
baryons or a caustic surface for the collisionless dark mat¬ 
ter (or galaxies). The position and velocity of the shock for 
a 7 = 5/3 ideal gas is very close to that of the outermost 
caustic surface for the dark matter; together they define a 
unique radius (commonly referred to as the virial radius 


Tvir) within which the mean enclosed density is ^ 100 times 
the background value. Since all matter outside the virial sur¬ 
face is infalling with the cosmic mix of components, then, by 
continuity, the baryon mass fraction measured at the virial 
radius must be unbiased fb{rvir) = 

The realistic case differs from the spherical one in sev¬ 
eral respects. Clusters forming in a fully three-dimensional, 
hierarchical fashion experience lumpy, asymmetric accretion 
directed by connecting filaments (West, Dekel & Oemler 
1987; Frenk et al. 1990; Evrard 1990a,b; Kang et al. 1994; 
Navarro, Frenk & White 1995; Tormen 1996). In addition, 
the ability of the intracluster gas to lose entropy via radia¬ 
tive cooling or increase it via feedback from galactic winds 
calls into question lessons learned assuming gravitationally 
induced shocks are the only entropy changing mechanism. 

These issues led White et al. (1993) to explore extreme 
models for dynamical baryon enhancement. They examined 
the evolution of infinitely dissipational (pressureless) gas 
and dark matter using both a spherical model based on 
Bertschinger’s (1985) solutions and three dimensional, gas 
dynamic simulations. Following their notation, define T((5c) 
as the ratio of enclosed baryon fraction within radius rsc to 
the cosmic value 

Mrsc) = r{Sc) (1) 

For the extreme case of a zero temperature gas, the White 
et al. three dimensional simulations show baryon enhance¬ 
ments smaller than the spherical case. For example, at 
Sc = 500, the spherical model predicts T(500) = 1.6 while the 
average of twenty simulations is T(500) = 1.25 and the maxi¬ 
mum simulation value is 1.5. Although the simulation results 
strictly apply to the case of a standard CDM initial fluctu¬ 
ation spectrum, results for other cosmologically reasonable 
Dm = 1 power spectra should not differ substantially, as the 
sensitivity of cluster dynamical histories to spectral shape is 
fairly mild (Lacey & Cole 1993; Crone, Evrard & Richstone 
1994; Navarro, Frenk & White 1996). 

When a realistic gas equation of state is used, com¬ 
bined N-body and gas dynamic simulations generally show 
the gas to be slightly more extended than the dark matter 
(Evrard 1990a; Thomas & Couchman 1992; Cen & Ostriker 
1993; Kang et al. 1994; Metzler & Evrard 1994; Navarro, 
Frenk & White 1995; Lubin et al. 1996), though the evi¬ 
dence is not universal (Anninos & Norman 1996). Energy 
transfer between these components during major mergers is 
a likely physical explanation (Navarro & White 1993; Pearce, 
Thomas & Couchman 1994), though the origin may be un¬ 
related to mergers (Chieze, Teyssier & Alimi 1996). The ex¬ 
tended gas structure implies a weakly rising baryon fraction 
with radius /gas(r) with rj^ 0.1 — 0.2 near the virial ra¬ 
dius, and a modest, overall baryon diminution T(500)^0.9 
(Frenk et al. 1996). The mild rise of gas fraction with radius 
appears consistent with the observations discussed below. 

2.2 Cluster mass estimation 

The bias and variance in the mass estimates inferred from 
observations have also been calibrated by numerical experi¬ 
ments of cluster formation incorporating gravity and gas dy¬ 
namics. A number of independent experiments over the years 
(Evrard 1990a,b; Tsai, Katz & Bertschinger 1994; Metzler & 



The Intracluster Gas Fraction in X-ray Clusters 3 


Evrard 1994; Navarro, Frenk & White 1995; Schindler 1996; 
Evrard, Metzler & Navarro 1996; Roetigger, Burns & Loken 
1996) show that, when exercised judiciously, accurate bind¬ 
ing mass estimates can be made with the standard, /3-model 
approach (Cavaliere & Fusco-Fumiano 1976). “Judicious ex¬ 
ercising” here means avoiding the cores of clusters where 
cooling flows and other complications occur (Tsai, Katz & 
Bertschinger 1994), avoiding clusters engaged in obvious ma¬ 
jor mergers (Roetigger et al. 1996), and avoiding extrapo¬ 
lating to very large radii where an equilibrium assumption 
is not justifiable. 

From an analysis of gas velocity moments, Evrard et 
al. (1996, hereafter EMN) propose rsoo as a conservative 
choice for the outer hydrostatic boundary of clusters. Using 
a sample of 56 cluster simulations evolved in different cos¬ 
mological backgrounds, and with a subset including input of 
mass and energy from galactic winds, they show that, within 
this radius, the gas is very close to hydrostatic, with a mean, 
mass weighted, radial Mach number of only a few percent 
(see Tables 3 and 4 of EMN). Binding mass estimates based 
on the standard, /3-model approach are nearly unbiased and 
have an intrinsic scatter of ~ 30% at rsoo. This compares 
well with the 15% rms deviation quoted by Schindler (1996) 
and the tens of percent variations seen in the experiments 
of Roettiger et al. (1996). 

In addition, EMN hnd that the variance in the mass 
estimates can be considerably reduced by eliminating the 
/3 parameter entirely [i.e., ignoring the X-ray image) and 
deriving the binding mass directly from the global, emis¬ 
sion weighted temperature Tx. The resulting scaling rela¬ 
tions are consistent with virial equilibrium expectations at 
a fixed density contrast T ~ GM/r ~ ScpcF. Calibrating 
the relation with the 56 experiments at Sc = 500 yields 

For the numerical sample, a standard deviation of only 15% 
in mass estimates results from this relation. 

Gas masses can be recovered with typically higher ac¬ 
curacy than the dark matter (Metzler & Evrard 1994; Roet¬ 
tiger et al. 1996), but this is based upon the assumption that 
the gas is smoothly distributed within intracluster space, not 
bound into filamentary or knotty clumps. There is no well 
defined physical model for such dumpiness, rather the mo¬ 
tivation for such a model comes from a desire to minimize 
the gas fraction in clusters while still providing the observed 
emission measure for bremsstrahlung. Since the latter scales 
as Pgas while the former as Pgas, a large, local ’’clumping 
factor” C =< p^as /pgas could reduce the amount of 

gas required to produce a fixed X-ray emission by a factor 
-C. 

There is some observational evidence arguing against 
such clumping. A very high signal-to-noise X-ray image 
of the Coma cluster shows no signs of fluctuations in the 
emission apart from that associated with individual galaxies 
(White, Briel & Henry 1993). In addition, measurements of 
the Sunyeav-Zel’dovich (SZ) decrement in clusters are con¬ 
sistent with expectations based on no significant clumping 


for reasonable values of the Hubble constant (Birkenshaw, 
Hughes & Arnaud 1991; Herbig et al. 1995; Carlstrom, Joy 
& Grego 1996). For example, again in Coma, Herbig et al. 
(1995) derive a Hubble constant Ho = 7lt25 km s“^ Mpc”*^ 
from OVRO observations of the SZ effect combined with 
the X-ray model of Hughes (1989) for the intracluster gas, 
which assumes no clumping. However, given the present un¬ 
certainty in these and other Hubble constant determinations 
(Kennicutt, Freedman & Mould 1995), there appears room 
for a systematic effect from clumping at the tens of percent 
level, and perhaps higher. This and other systematic effects 
are discussed in §4.3 below. 


3 THE MEAN CLUSTER GAS FRACTION 

In this section, I determine the mean gas fraction within 
rsooiTx) — calibrated by the numerical experiments, 
eq’n(||) — for the observational samples presented by White 
& Fabian (1995) and David, Jones & Forman (1995). 

White & Fabian (1995; hereafter WF) present gas frac¬ 
tions for 19 clusters derived from archival Einstein IPG ob¬ 
servations and temperatures compiled from the literature, 
many from David et al. (1993). They quote values for fgas 
at both a fixed metric radius of 0.5 h~^ Mpc and the cluster 
X-ray radius rx which is governed by the image quality. To 
evaluate gas fractions at rsoo{Tx), I use a mild, power law 
extrapolation of the data quoted at rx 

fgasirsOoiTx)) = fgas{rx) ( 4 ) 

with t/ = 0.17. This value of p is derived from the numerical 
simulations of EMN, and is consistent with the rise of gas 
fraction with radius in both the WF and DJF samples. The 
mean values of p in the WF data, derived by comparing 
the baryon fractions at 0.5Mpc and rx, is 0.13. This 
is biased somewhat low by the few clusters with short lever 
arm — removing 3 clusters with rx < 0.6 h~^ Mpc yields a 
mean 7/ = 0.15, with all but one value in the range 0.05 — 0.28. 

The resultant gas fractions at r 5 oo(Tx) for the WF sam¬ 
ple are plotted against cluster temperature as Hlled circles 
in Figure |^. Error bars on the gas fraction are 90% conh- 
dence limits and include contributions from the binding and 
gas mass errors. The errors quoted in Table 2 of WF in¬ 
clude only the gas contribution. To estimate binding mass 
uncertainty, I assume a fractional error equal to that in clus¬ 
ter temperature. Because the latter is generally asymmetric, 
so also is the binding mass error. Temperature values and 
errors are taken from the deprojected values in Table 1 of 
WF, except for two cases of very hot clusters — A2142 and 
A2163 — which have unusually large allowed lower temper¬ 
ature ranges from the deprojection method. For these, I use 
the “reference” lower bound on temperature quoted by WF 
rather than the deprojected values. 

The temperature errors are comparable in magnitude 
to those of the gas; average Ict fractional errors are 8.5% 
for the gas and 6.6/12.8% for the upper/lower temperature 
uncertainties. Because the gas emissivity in the IPG band 
is fairly insensitive to temperature over most of the range 
spanned by the WF data, the derived gas masses are nearly 
independent of temperature. Errors in gas and total masses 
can then be assumed uncorrelated, implying the fractional 


rsoolTx) = (1.24 ±0.09) 

MsoolTx) = (1.11T0.16) 
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Figure 1. The intracluster gas fraction within rsooCTx) from 
eq’n(^ plotted against X—ray temperature for the observational 
samples of WF (solid dots) and DJF (crosses). Errors in the tem¬ 
perature, omitted for clarity, are fractionally comparable to those 
in gas fraction. Estimates of the samples’ cumulative parent prob¬ 
ability distributions, derived from the normalized sum of the indi¬ 
vidual cluster asymmetric Gaussian contributions, are shown by 
the solid (WF) and dotted (DJF) lines. 

square error in fgas is the sum of the squares of the fractional 
errors in the gas and total masses. 

David et al. (1995; hereafter DFJ) provide in their Fig¬ 
ure 5 the gas fraction as a function of overdensity Sc directly. 
They present data on a range of systems, from elliptical 
galaxies to clusters; I use only the seven groups and clus¬ 
ters with Tx > 1 keV; namely, A539, A262, A2589, A2063, 
A1795, A85 and A2029. Errors in the binding masses for 
these systems are again taken to be proportional to the tem¬ 
perature errors listed in Table 1 of DJF. Gas mass errors are 
not quoted in their Table, for these I assume a fixed Icr error 
of 8.5%, the mean of the WF sample. This is likely to be an 
overestimate of the actual typical error, given the improved 
image quality of the ROSAT PSPC over the Einstein IPC. 
The results are insensitive to the choice of assumed error 
over the range 0 — 15%. The data are plotted as crosses in 
Figure ^ 

The cooler ROSAT sample appears to be more gas rich 
than its hotter IPC counterpart. Non-parametric statistical 
tests of location, such as the Mann-Whitney test and the 
Wilcoxon Signed-rank test, indicate sample inconsistency at 
better than 99% confidence for the raw data values. A simple 
way to incorporate errors is through estimates of the sam¬ 
ples’ parent probability distributions, derived from an un¬ 
weighted sum of the individual cluster asymmetric Gaussian 
contributions. The cumulative version of these are shown for 
each sample as the solid and dotted lines in Figure 

The range of likely baryon fraction values, as measured 
by the 5 to 95% confidence intervals, is large, 4 — 10 
and the two data sets provide consistent estimates for this 
range. However, they differ substantially in their median 
values of 5.7(WF) and 7.7 (djf) percent, re- 
fleeting the fact that the distributions are asymmetric and 
oppositely skewed. A likelihood analysis of the data yields 


(%, h-0.5) 
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Figure 2. Probability distribution function of the sample bi¬ 
weight mean gas fraction at rsoo obtained from bootstrap estima¬ 
tion. The left panel employs the original binding mass estimates 
from the literature while the right panel replaces the original mass 
estimates with values derived from and the virial scaling re¬ 
lation, eq’n(^. The dashed lines show the 95% confidence region 
for /gas(r50o) derived from the combined, revised samples. 

most likely values of 5.6 and 7.3respectively, 

very similar to the median values. 

A more robust assessment of sample location is made 
by employing a bootstrap procedure to estimate the range 
of sample means consistent with each data set. A bi-weight 
estimator of location is used, because of its superior per¬ 
formance for small data samples (Beers, Flynn & Gebhardt 
1990), but results using an ordinary mean or median are sim¬ 
ilar. The procedure creates a large number of trial samples 
with replacement, assuming data values distributed in an 
asymmetric Gaussian fashion. For each trial, the bi-weight 
mean is estimated and the resulting distribution over many 
trials constructed. 

The results, shown in panel (a) of Figure |^, are very 
nearly symmetric and Gaussian, the DJF sample less so be¬ 
cause of the influence of A2063, an outlier with fgaairsoo) = 
4.9 (see Figure |^. The bi-weight means and stan¬ 
dard errors of the samples — 5.82 ± 0.45 (WF) and 

7.77 ± 0.69h“^'^^% (DJF) — differ at the 2.4cr level, where 
the significance quoted is \xi — X 2 \l \/appropriate 
for independent Gaussian distributions with means ®i and 
X 2 and standard deviations ai and a 2 . This is supported by 
the appearance of the distributions in Figure |^( a). 

The direction of the discrepancy — lower temperature 
DJF groups having higher gas fractions than the richer clus¬ 
ters of the WF sample — is surprising, because adding the 
galaxy contribution to the baryon mass would serve to ex¬ 
aggerate the difference, not correct it. Adding the galax¬ 
ies would imply significantly larger baryon fractions in poor 
groups compared to rich clusters within the virial radius, an 
outcome not anticipated in hierarchical clustering scenarios. 
In fact, the opposite is the more likely expectation. Feed¬ 
back from galactic winds drives baryons preferentially out 
of low temperature systems (Yahil & Ostriker 1973; White 
1991; Metzler & Evrard 1994). It is possible that the discrep¬ 
ancy arises from sample selection criteria, but the fact that 
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neither sample is statistically well defined makes it difficult 
to address this issue. Qualitatively, it is difficult to under¬ 
stand why the generally lower quality IPC imaging should 
be biased in favor of lower gas fraction clusters. 

A possibility addressed here is that modest biases in the 
binding masses are responsible for the discrepancy. Though 
both are based on an underlying assumption of equilibrium, 
WF and DJF use slightly different approaches for estimat¬ 
ing binding masses. Only a small (~ 25%) systematic effect 
is required to bring the two sample means into agreement. 
Analysis of the numerical sample of EMN indicated the virial 
scaling relation provided a more accurate — in the sense of 
minimizing variance — mass estimator than the standard, 
/3-model method. If real galaxy clusters satisfy virial equi¬ 
librium to the deCTee established in the simulations, then ap¬ 
plication of eq’n (m to the observations should produce sim¬ 
ilarly accurate mass estimates. This suggests a re-analysis of 
the data, employing revised binding mass estimates derived 
from eq’n and the measured cluster temperatures T^. 

The result of repeating the bootstrap procedure with 
the revised binding masses is shown in Figure |^b). Appro¬ 
priate corrections have been made to maintain estimates at 
a density contrast 5c = 500. The WF data shift by a small 
amount, indicating good agreement between the virial scal¬ 
ing relation and the original deprojection mass estimates, 
and the variance decreases slightly. The DJF data shift more 
substantially to lower gas fractions. The direction and mag¬ 
nitude of the shift can be traced directly to the mean value of 
(3, which enters linearly in the original mass estimates. The 
seven groups and clusters of the DJF sample have /3 = 0.63, 
about 25% lower than the “magic” value of 0.79 inferred 
from inserting the radial scaling, eq’n (^), into the /3-model 
estimator (see EMN). 

The revised distributions provide consistent estimates 
of the mean cluster gas fraction between the two samples. 
Combining the data sets results in a highly statistically ac¬ 
curate estimate of the bi-weight mean gas fraction within 
’’soo for X-ray clusters 

fgasirsoo) = (0.060 ± 0.003) (5) 

where the quoted error is Icr. For comparison, the same 
limits derived from combining the original sample data is 
0.063 ± 0.004. In aligning the two samples, the revised mass 
estimates provide a small reduction in the uncertainty while 
not significantly affecting the location of the mean. The anal¬ 
ysis which follows employs the revised gas fraction value, 
but clearly similar numbers result if the original data are 
employed. 

This analysis of the population mean should not be mis¬ 
interpreted as presenting the value for the gas fraction within 
rsoo for all clusters. Simulations indicate intrinsic variations 
in the gas fraction at about the 15% level arise naturally 
from the different dynamical histories/states of a coeval pop¬ 
ulation (White et al. 1993; Cen & Ostriker 1993; Kang et 
al. 1994). Galactic feedback can cause gas loss by subsonic 
winds after cluster formation or by hindering collapse itself 
through early preheating. Though nearly all the 26 clusters 
in Figure have 90% confidence limits which overlap the 
limits in eq’n(^), there are clusters which have significantly 
less gas. An example is A576, which has a gas fraction lower 
by about a factor 2 than the mean (Mohr et al. 1996). Com¬ 
pact groups are extreme in this regard (Ponman et al. 1996; 


Pildis, Evrard & Bregman 1996). Loewenstein & Mushotzky 
(1996) present two cool clusters with different baryon frac¬ 
tions within their virial regions. Intrinsic variations in clus¬ 
ter gas content are thus both expected and observed, with 
the caveat that large amplitude variations go in only one 
direction, that of reducing gas content. There is currently 
no empirical evidence or theoretical justification supporting 
large baryon enhancements near the virial radius in clusters. 

4 IMPLICATIONS FOR Qm 

The arguments presented in §2 indicate that the mean 
baryon fraction within Sc = 500 is enhanced by at most a 
factor 1.25 by dynamical means during hierarchical cluster¬ 
ing and is more likely slightly below the universal value. 
This leads to two avenues for using the mean intracluster 
gas fraction to constrain Dm which are considered here. 

4.1 Upper limit approach 

Because the contribution of baryon sources add linearly, the 
mean baryon content of clusters must exceed the mean gas 
content fb > fgas. The mean gas fraction can thus be used 
to place an upper limit on Dm/Dj in eq’n(^) 

^ < T(4) f^cMrSc). (6) 

From the perspective of maximizing Dm, it is appropriate 
to take the 5%-ile lower limit on fgas along with the largest 
possible baryon enhancement. The simulations of White et 
al. yield an average enhancement T(500) = 1.25, and it is 
reasonable to assume this is appropriate for the effect on 
the population mean. The upper bound on the ratio of total, 
clustered to baryonic mass density is then 

^ < 23.1 (7) 

lib 

Comparison between the light elemental composition 
of the universe and primordial nucleosynthesis expectations 
places an upper limit on the mean baryon density. This can 
be combined with the above to express an upper limit on 
Dm directly. The exact value of the baryon density derived in 
this way continues to be a subject of current debate (Krauss 
& Kernan 1994; Copi, Schramm & Turner 1995; Steigman 
1995; Sasselov & Goldwirth 1995; Hata et al. 1996; Tytler, 
Fan & Buries 1996; Rugers & Hogan 1996). However, a con¬ 
sensus view is that a firm upper limit Dj < 0.026 h~^ exists 
simply from the abundanace of fragile deuterium in the hos¬ 
tile environment of the local solar neighborhood (Linsky et 
al. 1995). This value produces the constraint 

< 0.60. (8) 

The Einstein-deSitter case Dm = 1 requires a very low Hub¬ 
ble constant /i<0.36 (Bartlett et al. 1995). 

4.2 Best estimate approach 

This upper limit is generous in that it : (i) ignores the contri¬ 
bution of galaxies to the baryon fraction; (ii) assumes max¬ 
imal, dynamical baryon enhancement through use of an un¬ 
physical gas equation of state and (iii) uses the largest real¬ 
istic value of D(, combined with the smallest allowed value of 
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fgaa{r 5 Qo). An alternate perspective is to use “realistic” pa¬ 
rameter values to provide a “best” estimate of For this, 
I will assume : (i) galaxies make a 20 contribution rel¬ 

ative to the gas fbirsoo) = (1 + 0.2 h^^'^)fgas{r 50 o) as appro¬ 
priate for Coma (White et al. 1993) and (ii) a realistic equa¬ 
tion of state for the gas leads to a modest baryon diminution 
T(50(b = 0.85. Utilizing the central value of /gas(rsoo) from 
eq’n (ffl results in 

^ = (11.8 ±0.7) (11.8 ± 0.7) (9) 

l + 0.2h^^^ 

which is a factor of two smaller than the generous upper limit 
in eq’n (^). The Icr error quoted above is derived from the 
propagated statistical error of fgaa- The odd-looking slope 
of 4/3 in the second expression is a power law approximation 
to the 1.2/i®^^/(l ± 0.2h®/^) term; it is accurate to 2 percent 
over the range h £ [0.45, 1]. 

The final step to estimate flm requires a value for Qb- 
Recent inferences of the primordial deuterium abundance 
from quasar absorption line spectra produce two different 
preferred values, a low value = 6.2±0.8 x 10“® (Rugers 
& Hogan 1996) and a high value = 0.024±0.006 (Tytler, 
Fan & Buries 1996) derived from their respective high and 
low estimates for D/H. These values result in estimates and 
la errors of 


= 0.07 ±0.01 

high D/H, 

(10) 

= 0.28 ±0.07 

low D/H. 

(11) 


The latter value should be preferred over the former for sev¬ 
eral reasons. First, a recent analysis of high resolution Keck 
spectra of Q0014±813 — the best “high D/H” candidate — 
by Tytler, Buries & Kirkman (1996) provides no support for 
a level of high deuterium absorption, and suggests hydrogen 
interlopers as a likely explanation for the data. In addition, 
the value f2m'~0.3 is currently concordant with large-scale 
structure observations {e.g., Ostriker & Steinhardt 1995), 
while the low value <0.1 deduced from the Rugers & 
Hogan D/H estimate is decidedly difficult in this regard. 

4.3 Systematic effects 

The modest statistical error in these estimates is deceptive, 
since there are systematic uncertainties in the steps involved 
in the derivation. In the “minimal” approach adopted above, 
there is perhaps a factor two uncertainty in the contribution 
of galaxies relative to gas, implying about a 20% error in the 
baryon fraction estimate. In addition, there is approximately 
15% uncertainty in the appropriate value of T(500). Conser¬ 
vatively adding these contributions leads to an estimate of 
the overall systematic uncertainty of 35%. In addition, there 
are systematic effects in the nucleosynthesis determination 
of Di, which would enlarge the quoted statistical error (Au- 
douze, Olive & Truran 1997), but values Qbh^ > 0.024 are 
unlikely because of additional constraints from '^He, ®Li and 
^Li abundances {e.g., Lemoine et al. 1997). 

The actual error in this analysis could be larger if other 
systematic effects play a significant role. Current possibili¬ 
ties include the following. 

(i) Multi-phase intracluster gas — One can imagine a 
multi-phase structure, with cooler gas perhaps entrained 
along loops of magnetic field surrounded by hotter, more 


dilute plasma. In order to be competitive with thermal pres¬ 
sure and create a signihcant mass estimate bias through 
clumping, the magnetic pressure must be close to the ther¬ 
mal pressure outside the loops and the mass filling factor 
within the loops should be large. No specific model for the 
origin and maintenance of such a configuration has been pro¬ 
posed. On energetic grounds alone, it is difficult to imagine 
galaxies supplying enough magnetic held, and observations 
of Faraday rotation in background sources indicate that any 
strong helds must have small coherence lengths (Kim et al. 
1990; Kronberg 1994). Another energetic argument against 
signihcant magnetic pressure is the fact that the specihc 
thermal energy of the intracluster gas is very nearly equal 
to the specihc kinetic energy of the galaxies (Jones & For¬ 
man 1984; Lubin & Bahcall 1993), as expected if the thermal 
and kinetic pressures respectively support each component 
within the same potential well (the assumption underlying 
the /3-model of Cavaliere & Fusco-Femiano 1976). Spatially 
resolved X-ray spectroscopy, particularly of line emission in 
cooler (2 — 6 keV) clusters, will provide tests of such models, 
which are broadly similar to multi-phase cooling how mod¬ 
els (Sarazin 1996). Limited information is already available 
in broad-beam colors (Henriksen & White 1996). Data from 
ASCA, SAX, and soon XMM and AXAF, coupled with real¬ 
istic, dynamical and thermodynamical modeling (Teyssier, 
Chieze & Alimi 1996) should place stringent constraints on 
such multiphase models. 

In addition, since the electron pressure structure in the gas 
in a multi-phase model will diher from the standard, un¬ 
clumped model, Sunyaev-Zel’dovich measurements can be 
used to provide independent constraints on clumping. Roet- 
tiger’s (1996) gas dynamic merger simulations show that er¬ 
rors in SZ Hubble constant determinations are small in the 
standard scenario, providing hopeful prospects of detecting 
a modest signal from dumpiness. Present data probably can¬ 
not rule out clumping factors as large as 50%, but factors 
> 2 seem unlikely. 

(ii) Mass estimate errors — The simulations may be giv¬ 
ing a misleadingly simple picture of mass estimate accuracy. 
The good agreement of independent codes employing diher- 
ent gas dynamic methods on the same problem points the 
hnger at missing physics, rather than numerical inaccuracy, 
if a signihcant effect is to be found. Neither feedback from 
galactic winds (EMN) nor the inclusion of non-equilibrium 
thermodynamics in the intracluster plasma (Teyssier, Chieze 
& Alimi 1996) signihcantly alter the virial scaling relation, 
eq’n (^) . The fact that this binding mass estimator is able to 
reconcile the modest WF and DJF sample difference can be 
taken as a measure of empirical support for the simulation 
results, but that could be a misleading interpretation. Inde¬ 
pendent measures of cluster mass, particularly weak gravita¬ 
tional lensing (Tyson, Valdes & Wenk 1990; Kaiser & Squires 
1993) are needed to provide additional measures of mass es¬ 
timate accuracy. Comparison of the two methods near the 
virial radius in a few clusters indicates consistency within 
modest (~ 50%) statistical error ranges (Squires et al. 1996). 
Extending such studies to larger, statistical samples is im¬ 
perative. 

(in) Extrapolation errors — The X-ray data of the WF 
sample do not extend to rsoo, making extrapolation nec¬ 
essary. The extrapolation in gas fraction is modest, about 
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1 on average. As a further check on the employed 

procedure, the data values at rx can be used to predict val¬ 
ues at 0.5 h~^ Mpc, and the underlying distribution function 
constructed at that radius. The resultant 5, 50 and 95% con¬ 
fidence limits using the extrapolated gas fraction of 9.7, 13.8, 
and 22.9%, respectively (quoted for h=0.5), agree well with 
the directly determined values of 10.0, 13.8 and 22.3 from 
WF (their Figure 6). Significant errors from extrapolation 
are thus very unlikely. 

(iv) Hot dark matter — A sea of massive, light neutrinos 
would cluster differently from the cold dark matter typi¬ 
cally assumed in the dynamical simulations; their high en¬ 
tropy would prevent them from clustering in small potential 
wells. Experiments using viable cold plus hot dark matter 
(CHDM) models show that this effect is negligible in rich 
clusters at radii close to rgoo (Kofman et al. 1996). 

(v) Additional baryon contributions — Sources of baryons 
in clusters besides gas and galaxies exist, in the form of dif¬ 
fuse intracluster stars (Uson, Boughn, & Kuhn 1991; The- 
uns & Warren 1996) and MACHOS (Gates, Gyuk, & Turner 
1995; Alcock et al. 1996). The latter have not been directly 
detected in intracluster environments, and they may be sim¬ 
ply connected to the former. Regardless, extra baryons only 
drive the constraints on Qrn to lower values. Their contribu¬ 
tion probably does not exceed that of galactic starlight. 

(vi) Initial baryon fraction inhomogeneities — An early 
universe model which generates pre-inflationary baryon-to- 
total mass fluctuations and manages to preserve them into 
the post-inflationary epoch would circumvent the causal¬ 
ity argument mentioned in §2. Such a model would need to 
naturally couple high baryon overdensity to high mass over¬ 
density in order to produce an overestimate of the cosmic 
baryon fraction in clusters. Such a correlation would also 
avoid producing baryon deficient clusters, a comforting sit¬ 
uation since no large population of deep, “empty” potential 
wells exists (although Bonnet, Mellier & Fort (1994) have a 
candidate in the field of CL0024-I-1654). Consistency with 

= 1 would presumably be a natural feature of such mod¬ 
els. 

The sources of systematic uncertainty above (except the 
last) added in quadrature allow room for perhaps 70% ad¬ 
ditional upward error in Q.m. However, the magnitude of 
these effects — particularly of gas clumping for which there 
are no specihc, dynamical models — are currently not well 
understood. Forcing all effects in the same direction could 
probably manage consistency with Q.m = 1 and, in this case, 
the observational tests cited above should start uncovering 
the effects responsible in the near future. 


5 SUMMARY 

The gas fraction in clusters of galaxies provides informa¬ 
tion on the cosmic baryon mass fraction. Dynamical biases 
of clusters’ baryon content can be minimized by measuring 
masses near the virial radius, where gas dynamic experi¬ 
ments show equilibrium is valid and segregation processes in¬ 
efficient. At rsoo, the radius where the mean interior density 
is a factor 500 times the critical value, the bi-weight mean 
gas fraction of the combined, revised WF and DJF X-ray 


cluster samples is 0.060±0.003h“®/^. This value, when com¬ 
bined with our current understanding of cluster formation 
history and limits on from primordial nucleosynthesis, 
strongly favors D,m ~ 0.3 and rules out the possibil¬ 
ity of Qrn = 1 with high statistical significance, unless the 
Hubble constant is very low h < 0.4. These conclusions re¬ 
inforce, at greater statistical signihcance, earlier work based 
on different methods and data than those used here {e.g., 
Henriksen & Mamon 1994; Steigman & Felten 1995; Lubin 
et al. 1996). 

How serious is the case against Qrn = 1? A die-hard 
Einstein-deSitter advocate could espouse all the elements of 
the analysis leading to the upper limit in eq’n (^ and claim 
only a small (30% for h = 0.7) systematic effect is missing. 
However, this approach accounts for neither the hot. X-ray 
emitting gas nor the galaxies in clusters, and one might well 
be suspicious of an argument which ignores these two prin¬ 
cipal, observable components. The best estimate approach 
leaves one short by a factor > 3. It is possible that this gap 
is plugged not by a single large effect, but by several effects 
acting in concert, a situation reminiscent of the so-called 
“/1-discrepancy” in clusters (Evrard 1990b; Lubin & Bah- 
call 1995). There remains the possibility that a new element 
of X-ray cluster physics is simply missing from the present 
picture. Observational constraints on known sources of sys¬ 
tematic error should be vigorously pursued, along with more 
sophisticated theoretical modeling of intracluster plasma dy¬ 
namics and thermodynamics. 

The arguments presented here are independent of the 
value of the cosmological constant. If one favors a spatially 
flat universe, as motivated by simple models of inflation, 
obtained through a non-zero cosmological constant A, then 
the limits on the clustered matter component predict values 
of the required vacuum energy density Qa = A/3 Ho^ = 
1 — Qm- For example, for h = 0.7, the upper limit on 
Qm requires Qa > 0.28 while the best estimate approach 
gives Ha = 0.64 ± 0.09. The latter is consistent with the 
limit Ha < 0.66 derived from gravitational lensing argu¬ 
ments (Kochanek 1996) while marginally inconsistent with 
the Ha < 0.51 inferred recently by Perlemutter et al. (1996) 
from the magnitude-redshiftt relation for Type-la super¬ 
novae. 
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